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摘 要 运用 ANSYS 参数 化 设计 语言 建立 具有 真实 气动 外 形 和 复合 材料 铺 层 的 高 精度 有 限 元 过 结构 模型 ， 并 且 在 
Matlab 平台 调用 、 耦 合 基 于 BEM 理论 建立 的 气动 模型 ， 从 而 构建 风力 机 叶片 气 弹 分 析 模 型 . 充分 考虑 叶片 实际 受 载 情 
况 , 对 100 kW 风力 机 叶片 进行 弯 扭 而 全 加 载 气 弹 分 析 , 求解 了 耦合 与 非 耦 合 情 形 下 风力 机 的 响应 、 功 率 系数 和 弯 矩 等 ， 
总 结 不 同 风速 下 气 弹 看 合作 用 对 于 风力 机 的 影响 ,并 对 风力 机 叶片 设 i 30 
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Aeroelastic Coupling Ana se Blades Based on High Precision 
Rieé Element Model 


ZHAO Ying!’? Cai-Cail3 QIN Zhi-Wenl3 YANG Kel3 
(1. Institute of me Thermophysics, Chinese Academy of Sciences, Beijing 100190; 
2. a e School of Chinese academy of sciences, Beijing 100190; 
3. wd Sw key laboratory, Chinese academy of sciences, Beijing 100190) 
Abstract High precision finite element (FE) shell model are 《外 by Ansys parametric design 
language(APDL), real aerodynamic profile and composite layuBgsaie included in the model, and based 
on the Matlab platform, the aeroelastic model of wind t bihe(WT) blade is established by coupling 
traditional BEM theory and FE model. Taking the in-gs 
sponse, power coefficient and bending moment in ssotpling and uncoupling case of 100 kW WT with 


. 
load cases into account, the aeroelastic re- 


marizing the influence of aeroelastic coupling on 


bend-twist coupling blade loads are analyzed *5 
blade under different wind speeds, then, me improvements for the design of wind turbine blade. 
Key words finite element model; ba element-momentum theory; wind turbine blades; bend— 


twist coupling; aeroelastic coupling 


0 引 言 


随 着 风力 发 电 技术 的 快速 发 展 ， 风 电机 组 叶片 
长 度 不 断 增加 , 在 气动 设计 模块 , 传统 的 叶片 刚性 结 
构 模 型 假设 已 不 再 适用 ， 要 求 我 们 对 于 叶片 的 变形 
及 振动 作出 相应 假设 及 研究 〖 引 . 对 此 ,国内 外 学 
者 (6 开展 了 大 量 有 关 叶片 气 弹 模型 问题 的 研究 。 
李 本 立 等 中 采用 模 态 法 建立 了 风力 机 叶片 的 非 线 性 
模型 ; 储 程 等 四 以 Marc 通用 有 限 元 计算 程序 为 平 
台 , 通过 调用 FORCDT 载荷 子 程序 建立 二 维 气 弹 模 
型 ; Jong-Won Lee 等 外 改进 了 片 条 理论 , 用 有 限 元 
软件 建立 粱 单元 叶片 结构 模型 进行 耦合 计算 ; 莫 文 
收 稿 日 期 :2017-01-15; 修订 日 期 ，2017-03-22 
基金 项 目 : 国家 自然 科学 基金 资助 项 目 (No.51606196) 


威 等 0 构建 了 基于 多 体 动力 学 结构 模型 和 改进 型 
BEM 理论 气动 模型 相 耦合 的 气 弹 分 析 模 型 . 这 些 研 
究 中 的 气 弹 模 型 在 结构 模块 有 较 大 差别 ， 分 别 有 模 
态 法 、 多 体 动力 学 方法 和 有 限 元 方法 建 模 , 但 以 上 
结构 模型 没有 充分 考虑 真实 叶片 复合 材料 结构 铺 层 
特性 , 与 真实 叶片 结构 相差 较 多 ， 

因此 ， 本 研究 依托 我 单位 实验 风力 机 的 真实 载 
荷 特性 及 叶片 复合 材料 铺 层 特性 建立 结构 模型 ， 提 
升 模型 精度 ， 并 耦合 基于 BEM 理论 建立 的 气动 模 
型 ， 计 算 分 析 不 同 风速 下 叶片 气 弹 效应 的 影响 ， 为 


作者 简介 : 赵 座 (1991- ), 男 , 硕士 ,主要 从 事 风力 机 叶片 气动 弹性 的 研究 . zhaoying@iet.cn 
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4 期 
后 期 更 精确 的 研究 及 叶片 设计 及 优化 打下 基础 ， 
1 叶片 气 弹 模型 的 建立 
1.1 叶片 气动 模型 

在 风力 机 空气 动力 学 研究 中 ,计算 气动 性 能 常 
见 的 方法 主要 有 三 种 0 下， 一 种 是 通过 动量 理论 
与 叶 素 理论 结合 来 求解 风 轮 位 置 的 诱导 速度 ， 即 
“ 叶 素 - 动量 理论 ” (Blade Element-Momentum The- 
ory，BEM 理论 ); 第 二 种 是 通过 Biot-Savart 定律 
求解 叶轮 尾 流 涡 系 在 叶片 位 置 的 诱导 速度 ， 即 涡流 
理论 ; 第 三 种 是 通过 数值 计算 方法 求解 Euler 方程 
或 Navier-Stokes 方程 的 计算 流体 力学 方法 (Compu- 


输入 风力 机 、 风 场 基本 参 
数 ， 整 理 惨 型 升 阻 力 系数 、 


输入 截面 数 n， 邻 夫 1， 初 始 
诱导 因子 a、a 为 零 


计算 入 流 角 、 攻 角 


计算 Cu 和 C， 进 一 步 计 
算 修正 的 诱导 因子 wm、 


计算 Co 和 Cr， 积 分 求解 叶 
片 表面 的 作用 力 与 力矩 


图 1 气动 模块 大 代 流 程 图 


Fig. 1 Flow chart Iterative process of aerodynamic module 
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tational Fluid Dynamics，CFD)。 其中， 对 风力 机 叶 
片 进 行 气动 性 能 计算 的 方法 中 , 较为 简便 的 是 BEM 
理论 . 因为 BEM 理论 思想 简洁 明了 , 尤其 在 编程 运 
用 中 十 分 方便 ， 因 此 ， 在 研究 工作 中 得 到 了 广泛 地 
运用 中里 。 本文 以 BEM 理论 为 基础 ,加 入 叶 尖 、 叶 
根 损失 因子 进行 修正 , 运用 MATLAB 编写 程序 , 建 
立 风 力 机 叶片 的 二 维 气动 模型 ， 计 算 其 相应 的 气动 
性 能 ， 具 体 流程 如 图 1， 其 中 tol 为 精度 要 求 ， 此 处 
取 10-6. 
1.2 叶片 结构 模型 
结构 研究 方法 主要 有 三 种 : 一 种 是 基于 有 限 元 
理论 建立 模型 ; 第 二 种 是 基于 模 态 分 析 法 构建 模型 ; 
第 三 种 县 依据 多 体 动力 学 建立 模型 [4. 其中, 模 态 
过 程 简便 , 计算 量 小 , 但 精度 较 低 ; 有 限 


分 析 j 
建 模 过 程 复杂 ， 计 算 量 大， 但 对 于 复杂 模型 
原 度 高 ， 求 解 精度 高 ， 因 此 适用 于 要 求 较 高 的 


求解 过 程 ; 多 体 动力 学 不 论 模型 构建 还 是 求解 精度 
问题 都 居于 两 者 中 间 . 

本 文 为 了 准确 描述 叶片 复杂 的 复合 材料 结构 特 
性 ， 充 分 利用 有 限 元 分 析 软 件 ANSYS, 通过 参数 化 
语言 建立 叶片 三 维 结构 这 单元 模型 ， 并 根据 真实 叶 
片 复 合 材料 铺 层 情 况 进 行 铺 层 设置 ， 使 模型 具有 较 
好 的 仿真 程度 和 计算 精度 . 模型 结构 与 叶 根 部 位 某 


一 单元 的 铺 层 2 所 示 ， 其 中 单元 铺 层 共 36 层 ， 
图 中 显示 CN 20 层 . 
N 
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图 2 结构 模型 与 单元 部 分 铺 层 图 


Fig. 2 Structure model and partial layers of element 


2 计算 实例 


2.1 数值 计算 过 程 
通过 以 上 理论 与 软件 的 运用 , 针对 某 100 kW 风 


广 上 inmn”Avi 人 人 (下 甘 日 工 | 
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通过 BEM 理论 计 
攻 角 ， 求 解 升 
阻力 及 力矩 系数 


计算 诱导 因子 ， 
反复 迭代 至 满足 


变形 量 导 入 到 
气动 模块 ， 计 


对 前 后 两 次 变 
形 量 的 差 值 4 


对 结构 模型 做 
静 力 分 析 ， 提 


精度 和 和 用 


积分 求解 叶片 截 
面 的 力 及 力矩 


在 结构 模型 上 
加 载 


图 3 叶片 气动 弹性 分 析 过 程 


%g 
ve Fig. 3 Aeroelastic analysis process of blade 


力 机 ， 编 写 相应 的 MATLAB 与 APDL 程序 , 分 别 
建立 了 风力 机 叶片 的 二 维 气动 模型 与 三 维 结构 模型 . 
通过 弱 耦 合 方式 ， 建 立 了 风力 机 叶片 的 气动 弹性 分 
析 模 型 &1025， 其 过 程 如 图 3 所 示 . 


在 三 合 过 程 中 ， 通过 将 叶片 在 气动 力作 用 下 \\\、 
产生 的 扭转 变形 和 灾 曲 变形 反馈 到 叶片 的 气动 楼 硬 


中 ( 即 在 原 有 BEM 理论 的 攻 角 计算 式 中 气 弹 
扭转 角 的 影响 )， 从 而 反应 气 弹 耦 合影 响 ， 如 过 1 
所 示 。 


a=0—0 (1) 


式 中 , 9 为 叶 素 安装 角 , 0 = 8 十 p 十 (其 中 yp 为 叶 素 
扭 角 ; 6 为 桨 距 角 ; 6 为 叶片 在 考虑 气 弹 耦合 作用 后 
形成 的 扭转 变形 角 ， 简 称 气 弹 扭 转角 ). 
2.2 模型 验证 
2.2.1 气动 模型 验证 

基于 BEM 理论 ， 编 写 MATLAB 程序 进行 迭 
代 ， 并 与 bladed 软件 计算 结果 、phase VI 实 验 结 
果 09 进行 比较 , 如 图 4 所 示 . 图 4 中 点 从 左 到 右 风 
速 分 别 为 : 25 m/s、20 m/s、15m/s、13m/s、10m/s、7 
m/s, 可 以 看 出 , 在 叶 尖 速 比 较 小 时 误差 较 小 ,， 可 以 
将 其 用 于 分 析 7 m/s 左右 及 更 高 风速 时 的 工 况 。 


0.45 


3 
叶 尖 速 比 
图 4 功率 系数 与 叶 尖 速 比 关 系 


Fig. 4 _ Power coefficient variation with tip speed ratio 


2.2.2 结构 模型 验证 
结构 模型 由 APDL 参数 化 语言 建立 ， 通 过 

动 模型 计算 出 气动 力 ， 对 结构 模型 进行 加 载 ， 从 而 
计算 风力 机 在 气动 力作 用 下 的 响应 . 为 了 检验 结构 
模型 计算 的 准确 性 ,将 数值 计算 结果 与 静 力 实验 结 
果 [7 进行 对 比 .结构 模型 和 更 力 实验 加 载 位 置 相 
同 , 位 于 叶片 展 向 4m 与 8 m 处 . 二 者 叶 尖 挠 度 对 
比 显示 , 数值 计算 值 为 1.330 m, 静 力 实验 值 为 1.369 
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m， 相 差 2.8%。 这 主要 是 由 于 复合 材料 特性 的 离散 
性 导致 有 限 元 模型 和 实际 材料 存在 差别 ， 对 于 实际 
叶片 设计 ， 这 在 允许 的 误差 范围 内 ， 可 以 用 于 后 续 
的 结构 分 析 ， 
2.3 数值 计算 结果 

图 5 是 不 同 风 速 下 挥舞 方向 变形 量 与 叶片 径 向 
位 置 关系 ， 由 曲线 的 斜率 变化 可 以 看 出 叶 中 部 位 弯 
曲 程 度 最 为 严重 , 弯曲 弧度 大 , 拉 伸 程度 高 , 这 对 于 
叶片 本 身材 料 属性 形成 较为 严峻 的 考验 。 从 局 部 放 
大 图 可 以 看 出 叶片 变形 量 随 风 速 增 大 而 增加 ,在 叶 
尖 位 置 , 4 m/s 风速 时 考虑 气 弹 与 不 考虑 气 弹 造成 变 
形 量 相对 误差 是 1.48%，8 m/s 时 是 4.23%, 10 m/s 
时 是 4.58%。 因此 , 风速 升 高 加 剧 了 叶片 气 弹 耦合 响 
应 ， 随 着 风速 升 高 ， 入 流 角 增 大 ， 攻 角 也 随 之 增 大 ， 


O05 
造成 叶片 变形 量 增 大 , 但 当 风速 继续 升 高 , 风力 机 失 人 
束 ， 叶片 变 桨 ,导致 叶片 攻 角 减 小, 使 得 在 高 风速 区 XK 
叶片 的 变形 量 及 相对 误差 反而 逐渐 减 小 图 5 电站 


m/s 风速 时 相对 误差 是 2.85%, 24 m/s 时 是 39%/ 


挥舞 方向 变形 /m 


Se 2 
| | -~ 太 .23 m/s 无 气 弹 = 今 “ pp 
934 m/s 无 气 琐 -二 ,24 m/s 气 弹 |. 疡 


径 向 位 置 /m 
图 5 挥舞 方向 变形 量 与 径 向 位 置 关 系 


Fig. 5 Flapwise deformation variation with radial position 


叶片 变形 导致 的 扭 角 变 化 如 图 6, 其 中 , 扭 角 变 
化 大 的 区 域 主要 集中 在 叶 中 到 叶 尖 部 位 ， 因 为 此 处 
辟 型 气动 特性 较 好 ， 相 对 厚度 较 小 ， 在 弯 扭 耦合 作 
用 下 , 会 产生 较 大 的 气 弹 扭转 角 . 此 外 , 在 变 浆 运 行 
前 , 扭转 角 随 着 风速 升 高 而 增 大 , 在 4m/s 风速 下 叶 
尖 位 置 的 气 弹 扭转 角 为 0.533°, 6 m/s 时 为 0.749°, 8 
m/s 时 为 1.213°, 10 m/s 时 为 1.452"。 而 变 浆 后 ， 入 
流 角 依旧 升 高 , 攻 角 则 大 幅 减 小 , 进而 导致 扭转 角 减 
小 其 至 反 向 扭转 ， 最 终 产生 负 的 气 弹 扭 转角 , 在 15 
m/s 风速 时 为 0.516。, 20m/s 时 为 0.073°, 25 m/s 时 
为 -0.227°。 内 部 区 域 图 显示 了 风力 机 不 发 生变 桨 动 
作 时 扭 角 与 风速 关系 ,图 6 中 扭转 角 随 着 风速 升 高 
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不 断 增 加 ,最终 稳定 在 1.6” 左右 ， 由 图 6 可 得 , 不 
考虑 变 桨 时 , 随 着 风速 升 高 , 叶片 扭转 逐渐 变 大 , 即 
气 弹 斐 合作 用 随 着 风速 升 高 在 加 强 ， 


1.50 rl.50 
一 * 一 4m/s 人 ? 
125 “*%. 6ms E . 
1.25 上 
-8 m/s 让 
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4 Gs “+ 20 m/s A 的 
六 0.75 -一 .25 m/s Ea iA 
你 0.25 ma 站， 
E 0.50 了 1 内 
a 总 2 


4 6 
径 向 位 置 /m 
图 6 攻 角 变化 与 径 向 位 置 关系 


Fig. 6 Changeable angle of attack variation with radial 


Position 


图 7 是 风 功率 系数 与 风速 关系 曲线 , 图 中 虚线 
表示 的 是 考虑 气 弹 耦合 作用 后 的 风 功率 曲线 ， 而 实 
线 是 不 考虑 气 弹 耦合 作用 的 风 功率 曲线 。 因 为 考虑 
气 弹 耦 合作 用 造 咸 的 叶片 扭转 变化 后 ， 攻 角 会 发 生 
联动 变化 ,法 向 力 与 切 向 力 系数 也 随 之 改变 ， 从 而 
改变 叶 站 绪 存 ， 导 致 风力 机 功率 减 小 , 由 于 叶片 本 
oo 
恒 竹 略 ， 尤 其 是 在 7~15 m/s 风速 区 域 ， 其 中 在 7 


. Na 9 my/s、11 m/s 和 14 m/s 风速 时 功率 系数 分 别 


相差 0.04%、0.36%、1.51% 和 1.45%。 而 在 低 风 速 区 
域 (<7 m/s), 由 于 风速 较 低 , 叶轮 转速 也 较 低 , 载荷 
较 小 , 气 弹 效 应 不 明显 ; 在 高 风速 区 域 (>15 m/s)， 
由 于 叶片 变 桨 ， 导 致 叶片 攻 角 急剧 下 降 ， 乃 至 在 叶 
尖 处 产生 负 攻 角 ,， 在 此 区 域内 攻 角 改变 引起 的 法 向 
力 与 切 向 力 系数 变化 不 明显 ， 因 而 叶片 气 弹 看 合 效 
应 不 显著 。 

图 8 是 摆 振 方向 弯 矩 随 径 向 位 置 变 化 曲线 , 图 中 
分 别 对 比 了 考虑 与 不 考虑 叶片 气 弹 看 合作 用 对 弯 先 
的 影响 . 图 8 显示 , 随 着 风速 升 高 , 控 振 弯 逢 在 考虑 
与 不 考虑 气 弹 看 合作 用 后 的 差 值 增 大 ， 当 风速 达到 
一 定 程度 时 , 风力 机 开始 变 桨 动作 , 攻 角 减 小 , 摆 振 
弯 矩 之 间 的 差 值 反 而 减 小 , 其 中 在 6 m/s、8 m/s、 10 
m/s、15 m/s 和 20 m/s 风速 时 弯 矩 的 相对 误差 分 别 
是 2.53%、2.75 名 、2.80%、2.27% 和 1.68%。 由 此 可 以 
得 出 ， 随 着 风速 逐渐 升 高 ， 气 弹 耦 合作 用 效应 先 增 
大 而 后 减 小 。 
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图 7 功率 系数 与 风速 关系 


Fig. 7 Power coefficient variation with wind speed 
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图 8 不 同 风速 下 摆 振 方向 弯 矩 随 径 向 位 置 变 化 


Fig. 8 Edgewise bending moment variation with radius under 


different wind speed 


2.4 气 弹 设计 改进 
气 弹 耦合 效应 造成 的 转 变化 会 导致 功 率 路 
考分 不同 风 巡 下 转角 的 变化 全 中 ， 全 人 
的 所 转角 变化 什 , 从 而 对 叶片 的 设计 扭 角 进 行 改进 ， 
使 得 风 功 率 系数 在 考虑 气 弹 谭 全 效应 后 有 所 提高 ， 
改进 前 、 后 的 功率 系数 对 比如 图 9， 图 中 改进 后 的 


变化 前 
--0- 扭 角 变化 后 
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图 9 功率 系数 与 风速 关系 


Fig. 9 Power coefficient variation with wind speed 
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功率 输出 在 7 m/s 风速 前 基本 不 变 , 在 7 m/s 风速 
后 有 所 提升 , 其 中 8 m/s、9 m/s、11.3 m/s 和 15 m/s 
风速 时 功率 系数 分 别提 高 0.02%、0.32%、1.45% 和 
2.13%。 通过 预 扭转 改进 设计 扭 角 能 提高 风力 机 功率 
输出 ， 可 为 后 续 考 虑 气 弹 耦 合 效应 进行 叶片 优化 设 
计 提 供 参 考 ， 


3 结 


结 论 


本 文 发 展 了 一 种 叶片 气 弹 看 合 模型 ， 该 模型 运 

用 ANSYS 建立 叶片 复合 材料 铺 层 的 结构 模型 ， 耦 

合 基于 BEM 理论 建立 的 气动 模型 ， 形 成 叶片 气 弹 

分 析 模 型 ， 并 对 100 kW 风力 机 组 叶片 进行 气动 弹 
ed ts 

弹 耦 合 的 影响 先 随 着 风速 的 升 高 而 增 大 ， 


7 
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， 变 桨 程度 不 断 提 高 ， 气 弹 斐 合 的 影响 将 逐渐 
减 小 。 
2) 通过 气动 弹性 结果 分 析 可 知 ,在 正常 运行 风 
速 区 域 , 叶片 气动 弹性 效应 会 造成 叶片 功率 输出 下 
降 , 导致 达 不 到 设计 功率 。 因此 , 在 叶片 气动 设计 时 
要 充分 考虑 气 弹 效应 影响 ， 
3) 分 析 不 同 风速 下 气 弹 效应 造成 的 叶片 扭 角 
这 化 次 出 了 GR 交 天 功率 的 种 施 ， 关联 得 一 条 


效果 ， X. 
Cr 
献 


参 者 .》 
人 风力 机 气动 弹性 与 动态 响应 计算 [D]. 南京 : 南京 
航空 航天 大 学 , 2011 

CHEN Jiahui. Calculation of Aeroelastics and Dynamic 
Responses of Wind Turbine [D]. Nanjing: Nanjing Uni- 
versity of Aeronautics and Astronautics, 2011 


[2] 信 伟 平 . 风力 机 旋转 叶片 动力 特性 及 响应 分 析 [D]. 汕头 : 汕 
头 大 学 , 2005 


XIN Weiping. Dynamic Characteristics and Response 
Analysis of Wind Turbine Rotor Blades [D]. Shantou: 
Shantou University, 2005 

王旭东 , 陈 进 , Shen Wenzhong, 等 . 基于 气动 弹性 理论 的 风 
力 机 叶片 厢 合 分 析 [可 . 太阳 能 学 报 , 2010, 31(1): 96-99 
WANG Xudong, CHEN Jin, SHEN Wenzhong, et al. 
Wind Turbine Blades Coupling Analysis Based on the 


[3] 


Aeroelastic Theory [J]. Acta Energiae Solaris Sinica, 2010, 
31(1): 96-99 

[4 钱 若 军 , 董 石 腾 ， 袁 行 飞 等 . 流 固 耦合 理论 研究 进展 [可 . 空 
间 结 构 , 2008, 14(1): 3-14 
QIAN Ruojun, DONG Shilin, YUAN Xingfei, et al. Re- 
search Progress of Fluid-Solid Coupling Theory [J]. Space 
Structures, 2008, 14(1): 3-14 

[5] 陈 大 伟 , 杨 国 伟 ， 静 气动 弹性 计算 方法 研究 [可 . 力学 学 报 ， 
2009, 41(4): 469-479 
CHEN Dawei, YANG Guowei. Research of Static Aero- 


2 TY 


201705.00917V1 


chinaXiv 


4 期 


ey 


[11 


[12 


] 


elasticity Calculation Method [J]. Acta Mechanica Sinica, 
2009, 41(4): 469-479 

孙 良 . 弹性 变形 对 机 喜气 动 特性 影响 的 研究 [D]. 哈尔滨 : 哈 
尔 滨 工程 大 学 , 2012 

SUN Liang. Research on the Effects of Elastic Deforma- 
tion on Aerofoil Aerodynamic Characteristics [D]. Harbin: 
Harbin Engineering University, 2012 

李 本 立 , 安 玉 华 . 风力 机 气动 弹性 稳定 性 的 研究 [可 . 太阳 能 
学 报 , 1996, 17(4): 314-320 

LI Benli, AN Yuhua. Wind Turbine Aeroelastic Stabil- 
ity Study [J]. Acta Energiae Solaris Sinica, 1996, 17(4): 
314-320 

傅 程 , 王 延 荣 . 风力 发 电机 叶片 气动 弹性 响应 分 析 [可 . 机 械 
设计 与 研究 , 2009, 25(1): 68-70 

FU Cheng, WANG Yanrong. Wind Turbine Blade Aeroe- 
lastic Response Analysis [J]. Machine Design and Re- 
search, 2009, 25(1): 68-70 

Lee J W, Lee J S, Han J H, et al. Aeroelastic Analysis 
of Wind Turbine Blades Based on Modified Strip Theory 
[J]. Journal of Wind Engineering and Industrial Aerody- 
namics, 2012, 110: 62-69 

MO Wenwei, LI Denyuan, WANG Xianneng, et al. 


Aer 
lastic Coupling Analysis of the Flexible Blade of 下 
Turbine [J]. Energy, 2015, 89: 1001-1009 


Burton T, Sharpe D, Jenkins N , et al. es 
1: 46-65 


Handbook [M]. England: Antony Rowe 
张波 成 , 万 志 强 , 杨 超 . 基于 气动 三 气动 弹性 响应 


分 析 加 . 工程 力学 , 2011, 28(7): 
ZHANG Bocheng, WAN Zhigfia ANG Chao. Static 


\ 


赵 鹰 等 : 基于 高 精度 有 限 元 模型 的 叶片 气 弹 耦 合 分 析 


[13] 


[13] 


[4 


EC 


和 


[16] 


ChinaXiv 合 作 期 刊 


753 


Aeroelastic Response Analysis Based on Aerodynam- 
ics Three-Dimensional Force [可 . Engineering Mechanics, 
2011, 28(7): 217-222 
贺 德 馨 , 陈坤 , 张 亮 亮 等 . 风 工程 与 工业 空气 动力 学 [M]. 北 
京 : 国防 工业 出 版 社 , 2006: 87-91 
HE Dexin, CHEN Kun, ZHANG Liangliang, et al. Wind 
Engineering and Industrial Aerodynamics [M]. Beijing: 
National Defence Industry Press, 2006: 87-91 
唐 迪 , 陆 志 良 , 郭 同 庆 . 大 型 水 平 轴 风 力 机 叶片 气动 弹性 计 
算 [可 . 应 用 数学 和 力学 , 2013, 34(10): 1091-1097 
TANG Di,LU Zhiliang,GUO Tongqing. Large Horizon- 
tal Axis Wind Turbine Blade Aeroelastic Calculation [J]. 
Applied Mathematics and Mechanics, 2013, 34(10): 1091-— 
1097 
Carrion M, Steijl R, Woodgate M, et al. Aeroelastic Anal- 
ysis of Wind Turbines Using a Tightly Coupled CFD-CSD 
Methgd [J]. Journal of Fluids and Structures, 2014, 50: 
3 

军 . 不 同 风 况 下 风力 机 尾 流 特性 研究 [D]. 北京 : 中 国 
学 院 工 程 热 物理 研究 所 , 2014 
WANG Shengjun. Study on Wake Characteristics of 
Wind Turbine Operating in Various Wind Conditions [D]. 
Beijing: Institute of Engineering Thermalphysics, Chinese 
Academy of Sciences, 2014 
CHEN Xiao, QIN Zhiwen, YANG Ke, et al. Numerical 
Analysis and Experimental Investigation of Wind Turbine 
Blades with Innovative Features: Structural Response and 
Characteristics [J]. Science China Technological Sciences, 
2014, 58(1): 1-8 


